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ABSTRACT 

The black hole (BH)-bulge correlations have greatly influenced the last decade of effort to understand galaxy 
evolution. Current knowledge of these correlations is limited predominantly to high BH masses {Mbh^ 10^ Mq) 
that can be measured using direct stellar, gas, and maser kinematics. These objects, however, do not represent the 
demographics of more typical L < L* galaxies. This study transcends prior limitations to probe BHs that are an 
order of magnitude lower in mass, using BH mass measurements derived from the dynamics of H2O megamasers 
in circumnuclear disks. The masers trace the Keplerian rotation of circumnuclear molecular disks starting at radii 
of a few tenths of a pc from the central BH. Modeling of the rotation curves, presented by Kuo et al. (2010), 
yields BH masses with exquisite precision. We present stellar velocity dispersion measurements for a sample 
of nine megamaser disk galaxies based on long-slit observations using the B&C spectrograph on the Dupont 
telescope and the DIS spectrograph on the 3.5 m telescope at Apache Point. We also perform bulge-to-disk 
decomposition of a subset of five of these galaxies with SDSS imaging. The maser galaxies as a group fall below 
the Mbh-ct* relation defined by elliptical galaxies. We show, now with very precise BH mass measurements, that 
the low-scatter power-law relation between Mbh and cr* seen in elliptical galaxies is not universal. The elliptical 
galaxy Mbh - cr* relation cannot be used to derive the BH mass function at low mass or the zeropoint for active 
BH masses. The processes (perhaps BH self-regulation or minor merging) that operate at higher mass have not 
effectively established an Mbh - <^* relation in this low-mass regime. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert — galaxies: bulges 



1 . DYNAMICAL BLACK HOLE MASSES 

The past two decades have seen substantial improvements in 
our understanding of the demographics of nuclear supermas- 
sive black holes (BHs), predominantly driven by the availabil- 
ity of BH mass measurements for nearby galaxies from dy- 
namical techniques. The advent of the Hubble Space Tele- 
scope has enabled a few dozen dynamical BH mass mea- 
surements, which reveal that supermassive BHs are ubiq- 
uitous cornponen ts of bulge-dominated galaxy centers (e.g., 
iKormendvl |2004|) . Furthermore, there appears to be a re- 
markably tight correlation between BH mass and bulge prop- 
erties, including velocity dispersion (the Mbh ~ f^* relation; 



Gebhardt et al.'200rt'Ferr arese & Merritll2000l:lTremaine et all 
2002 : Giiltekin et al. 2009) , lum inosity (the MB H-^buige rela- 



tion, e .g., IMarconi & H unt 2003), and mass (e.g.. Harin g & Rixl 
|2004|) . The apparently tiny scatter in the Mbh-o"* relation has 
led to a common view that the growth of BHs is intimately con- 
nected to the growth of the surrounding galaxy (see discussion 
in[Ho 2004). In this paper, we explore the demographics of a 
new sample of dynamical BH masses measured using observa- 
tions of circumnuclear masers (Kuo, C. Y. et al. in prep.). 

Apart from our own Galactic Center, where indi vidual stars 
are used as test particles to weigh the BH (e.g., iGhez et al] 



|2b08; Gill essen et al.l 120091) . the most precise BH masses are 
derived using wat er megamaser spots in a Keplerian circum- 
nuclear disk (e.g.. lMivoshi et al]ll995b . Such megamasers are 
very close to an ideal dynamical tracer of the central mass at 
the galaxy center, as they delineate a disk with a well-measured 
Keplerian rotation curve within fractions of a pc from the cen- 
ter. They are detectable when the disk is aligned within a few 
degrees of the line of sight. Masers detected in edge-on nuclear 
disks have characteristic triple-peaked spectral profiles, mak- 
ing th e disk inclination angle well-known in these cases dLol 
l2005h . BH mass measurements based on the Keplerian rota- 
tion curves of circumnuclear megamaser disks are not prone 
to the same systematic uncertainties as the dynamical meth- 
ods based on optical observations. In the case of gas dynam- 
ical methods, it can be very difficult to determine the level of 
non-virial gas support (e.g., due to turbulence) but neglecting 
these extra sourc es of sup port can lead to significant bias in the 
final result (e.g., i Barth et al. 2001). Stellar dynamical meth- 
ods are sensitive to assumptions about galaxy anisotropy and 
radial variation in dark matter fraction t hat are not well con- 
strained observationally in general (e.g ., iGebhardt & Thornasl 
120091: Ivan den Bosch & de Zeeuw 2010). 

The cleanest and best-studied megam aser galaxy is the 
nearby Seyfert galaxy NGC 4258 (e.g., iNakai et al.1 1 1 9931: 
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Table 1. Sample and Observations 



Galaxy 


RA 


Dec. 


D 


Hubble Type 


Mb 


B-V 


Instr. 


Obs. Date 


PA 


^exp 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


NGC1194 


03:03:49.1 


-01:06:13 


52.0 


SAO+ 


-20.3 


1.03 


B&C 


2010-03-12 


139 


1500 


NGC2273 


06:50:08.6 


+60:50:45 


26.0 


SB(r)a: 


-20.2 


0.94 


DIS 


2009-05-01 


68 


1500 


UGC3789 


07:19:30.9 


+59:21:18 


50.0 


(R)SA(r)ab 


-20.5 




DIS 


2009-05-01 


165 


2700 


NGC2960 


09:40:36.4 


+03:34:37 


71.0 


Sa? 


-20.8 




B&C 


2010-03-12 


40,130 


3600,3600 


IC2560 


10:16:18.7 


-33:33:50 


41.8 


(R')SB(r)b 


-21.3 


0.97 


B&C 


2010-03-14 


44,134 


4800,4800 


NGC3393 


10:48:23.4 


-25:09:43 


53.6 


(R')SB(rs) 


-20.9 


0.88 


B&C 


2010-03-12 


158,68 


5400,4800 


NGC4388 


12:25:46.7 


+12:39:44 


19.0 


SA(s)b: 


-21.9 




B&C 


2010-03-13 


91,1 


3600,3600 


NGC6264 


16:57:16.1 


+27:50:59 


136.0 


S? 


-21.1 




B&C 


2010-03-13 


24,116 


3000,3600 


NGC6323 


17:13:18.1 


+43:46:57 


105.0 


Sab 






DIS 


2010-03-13 


176 


2100 



Note. — Col. (1): Name. Col. (2): Right Ascension (HH:MM:SS; J2000). Col. (3): Declination (DD:MM:SS; J2000). Col. (4): 
Distance (Mpc). Col. (5): Hubble type from RC3 (de Vaucouleurs et al. 1992) as compiled in NED. Col. (6): Absolute B-band magnitude 
(mag) from Hyperleda. Col. (7): B-V (mag). Col. (8): Instrument used for observations. The B&C spectrograph is on the Dupont 
telescope at Las Campanas, while DIS is on the 3.5m telescope at Apache Point Observatory. Col. (9): Date of observation. Col. (10): 
Position angle (°) East of North. When two slit positions were used, both PAs are listed. Col. (11): Total duration of spectroscopic 
exposure (sec). Two values correspond to major and minor axis exposure times respectively. 



GreenhiU et all 119951: iMivoshi et all 119951: iHerrnstein et"ai] 
19961 and for a complete literature review see Lo 2005). Obser- 



vations of centripetal acceleration from the frequency change 
of the maser components at the systemic velocity of the sys- 
tem, combined with VLBA images of the megamaser disk, al- 
low for the most p recise extragalactic an gular diameter dis- 
tance measurement JHerrnstein et al.ll999h . Other observations 
provide detailed constraints on the thickness, velocity struc- 
ture, and tor oidal magnetic field in an accretion disk (e.g. , 
Modiaz et a l] 120051: IHerrnstein et al.1 120051: lArgonet alj|2007t 
Humphrevs et al.l |2008[) . For a long time it proved difficult 



to de tect additional well-defined Keplerian m egamas er disks, 
(e^..Green hill et alJll997l:lBraatz etalJI 19971: iKondratko et all 
l2006b.) . but the situation has changed dramatically with the 
advent of the Green Bank Telescope (GBT) of the NRACB 
Thanks to the increased sensitivity of the GBT with a wide- 
band spectrometer, more megamaser candidates have been dis- 
covered (|BrMz&^Gu£liucdl2Q0M|&eeiMLeLay|2009|). Of 
the ~ 130 extragalactic megamasers, approximately half were 
discovered with the GBT and ^ 20 show evidence of originat- 
ing in a disk. 

The primary goal of ongoing megamaser searches is to find 
additional disk masers that can be used for angular diameter 
distanc e measurements i n the Megamaser Cosmology Project 
(MCP; iReid et al.ll2009t iBraatz et all [20101) . With a sufficient 
number of angular diameter distances to galaxies in the Hub- 
ble flow, the MCP aims to determine a pre cise independent 
Hubble constant t o a few percent accuracy dReidetal .1120091: 
iBraatz et al.ll2010 l). An important outcome of the MCP is that 
the Keplerian megamaser disks also provide a precise measure- 
ment of the enclosed mass within a few tenths of a pc from 
the center of these spiral galaxies, dominated by a BH at the 
center, a valuable resource in its own right. The BH masses 
for seven maser galaxies are presented in Kuo et al. (in prepa- 
ration). The present paper examines the demographics of the 
new megamaser galaxies, focusing on the shape and scatter of 
BH-bulge relations for these galaxies. 



To date, the vast majority of dynamical BH masses have been 
obtained in massive elliptical and SO galaxies (Fig. 1). There 
are only five targets with BH mass measurements < 10'' Mq 
(one of those being the Milky Way) and there are only 1 1 spi- 
ral galaxies with Hubble type of Sa or later in the sample of 
67 galax ies used in the most r ecent calibration of the Mbh - cr* 
relation JGiiltekin et al.l 120091 Fig. 1). The disproportionate 
representation of massive galaxies is easily understood. Local 
elliptical galaxies are relatively dust-free, and the BHs are mas- 
sive enough that we can resolve their gravitational spheres of 
influence at distances of tens of Mpc. Spiral galaxies tend to be 
dusty and star-forming, complicating dynamical measurements. 
Our knowledge of BH demographics at low masses is rudimen- 
tary, such that the shape and scatter in BH-bulge scaling rela- 
tions as well as the B H occupation fractio n in spiral galaxies 
remain uncertain (e.g.. lGreene & Holl2007h . The main focus of 
this work is to exploit the new sample of reliable BH masses 
afforded by the new maser studies to explore BH demographics 
for BHs with Mbh'^ 10^ Mq in spiral galaxies. 

The paper starts with a description of the sample (§2). We 
then describe the observations and data reduction (§3) and the 
spectroscopic (§4) and imaging (§5) analysis. In §6 we present 
the primary result of the paper, namely the location of the maser 
galaxies in the Mbh-c* plane, while in §7 the scalings between 
Mbh and bulge luminosity/mass are briefly discussed. The im- 
pact of our results on our knowledge of BH demographics, par- 
ticularly as probed with active galaxies, is discussed in §8, and 
we summarize and conclude in §9. 

Throughout we assume the following cosmological parame- 
ters to calculate distances: Hq = 70 km s"' Mpc"', f^m = 0.30, 
and f^A = 0.70. 

2. THE SAMPLE AND THE BLACK HOLE MASSES 



Over the years there have 
searches for m egamaser galaxies (e.g.. 



Inc. 



bee n a numb er of 

^_,^^__^^^^^____„ IBraatz etalJ 119971: 

IKondratko et"arj I2006bt Ifiiraatz & Gugliuccil 12008^ By 
and large, the surveys have targeted known obscured ac- 
tive galaxies, with a detection rate of ^5%. Once 

' The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, 
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Fi g. 1 . — The Mbh ""■* relation as defined by nearby inactive galaxies from 
IGUl tekin et al. (2009), color-coded by morphological type of the host galaxy. 
We distinguish between elliptical (red circles), SO (green triangles), and spi- 
ral (blue squares) galaxies. The distribution of stellar velocity dispersion is 
shown in the black histogram in the top panel, also divided into the contribu- 
tion of elliptical (red long-dashed lines), SO (green short-dashed lines), and 
spiral (blue dot-dashed lines) galaxies. Similarly, the Mbh distributions are 
shown in the right-hand panel. The histograms, as well as the black dashed 
lines, highlight the dearth of objects at it*< 125 km s"' and M^n< lO' Mq, 
the regime investigated by the current work. 



detected in a single-dish survey, megamasers that show disk 
characteristics, e.g., both systemic and high-velocity Doppler 
components, and that are bright enough, are subsequently im- 
aged with Very Long Baseline Interferometry (VLB!) to map 
the positions and velocities of individual maser spots. The 
map determines whether a maser is in a clean dynamical sys- 
tem or has more complicated structures such as outflows (e.g., 
iKondratko eTsHlQQS) . 

In this paper we present velocity dispersions both for 
the Kuo et al. samp le and two other objects, NGC 3393 
(IKondratko et al.ll2008l) and IC 2560 (Ishiharaet al. 2001). In 
the latter case, while it is assumed that the systemic masers arise 
from the inner edge of the masing disk, without VLBI mapping 
the enclosed mass is uncertain by factors of a few. We expect 
better data will yield a concrete mass. The VLBI (and thus 
mass constraints) for the remaining seven objects are presented 
for the first time in C-Y. Kuo et al. in preparation. The details of 
the BH mass fitting are described in Kuo et al. but are reviewed 
here briefly for completeness. The BH masses are derived by 
fitting a Keplerian rotation curve to the positions and veloci- 
ties of the maser spots. The distance is fixed using an assumed 
value of Hq and the inclination of the maser disk is assumed to 
be 90°. All of the megamaser disks for which we have mea- 
sured an inclination are within 6° of being perfectly edge-on, 
so the inclination contributes less than a 1 % error to the mea- 
surement of the black hole mass. In NGC 4388 we could not 
measure a disk inclination, but even if the disk is 20° from edge- 
on, the contribution to the error in the BH mass would be only 
12%, which is comparable to the BH mass error caused by the 
distance uncertainty (11%). Ongoing monitoring of the most 
promising targets will allow the measurement of accelerations 
in the systemic features. Joint modeling of Mbh and distance 



for these targets will l ead to a tighter constraint on each (e.g., 
iHerrnstein e"tanil999l) . 

It is interesting to note that the BH masses are strongly clus- 
tered around ^10^ Mq. We suspect that this is a simple con- 
sequence of picking random galaxies from the local active- 
galaxy mass function. Without accounting for incompleteness, 
the observed acti ve mass functio n is strongly peaked at a mass 
~ 10^ Mq (e.g., Heckman et al. 2004). At higher mass, there 
are very few systems radiating at an appreciable fraction of their 
Eddington limit and we are woefully incomplet e at lower mass 
dOreene & Holl2007HSchuize & Wisotzkill2010l) . While we do 
not sample the full range of BH masses in the local Universe 
with this sample, our selection is unbiased with respect to bulge 
mass and velocity dispersion and thus does not invalidate our 
results. 

The megamaser galaxies are found predominantly in early- 
to-mid-type spiral galaxies, ranging in Hubble type from SO 
(e.g., NGC 1194) to Sb (e.g., NGC 4388). Table 1 incl udes 
the RC3 morphological types jde Vaucouleurs et a l.lll992l ) and 
.B-ban d luminosities and colors from Hyperleda (Paturel e t alj 
l2003h . As noted above, the parent sample consists predom- 
inantly of known active galaxies, ranging in distance from 
18 to 150 Mpc. Only UGC 3789 did not have a docu- 
mented active galactic nucleus (AGN) when the circumnu- 
clear maser was discovered, although is it a narrow-line AGN. 
With the exception of N GC 4388, which has a weak broad 
Ha line (iHo et al.lll997h . the active galaxies are classified as 
narrow-line (obscured) AGNs from the optical spectra. Fur- 
thermore, those with X-ray spectra have larg e column densi- 
ties. Most are Compton thick (Nh > iO^ cm~^; Guainazzi et al 



'2005; Madeiski et al. 2006; Zhang et al. 2006; K ondratko et al 
.2006a; .GreenhiU et al...2008; .Zhang et al...2010.) . The only ex- 
ceptions are NGC 4388, which has been observed to undergo 
large variations in column density over time scales of hours 
(Elvis et al. 2004), and NGC 6264, NGC 6323, and UGC 3789, 
which have not yet been observed in the X-rays. The attempted 
XMM-Newton observations of the former two were destroyed 
by flaring. 

3. OBSERVATIONS AND DATA REDUCTION 

The primary observations discussed here were obtained us- 
ing the Dupont 2.5m telescope at Las Campanas Observatory 
in the South, and the Apache Point Observatory (APO) 3.5m 
telescope in the North. In the Appendix we describe an aux- 
iliary data set that we use to cross-check our velocity disper- 
sion measurements. At the Dupont telescope, we used the 
Boiler & Chivens (B&C) spectrograph with the 600 lines mm"' 
grating and the 2" slit, which yielded an instrumental resolu- 
tion of (Tr ^ 120 km s"' and a wavelength coverage of 3400- 
6600 A, covering the 4000A break, the G-band at 4304A, 
H/3h-[0 III] AA4959,5007, the Mg \b A A5 167, 5 173, 5 184 
triplet, and (in most cases) Ha. We observed at two slit po- 
sitions (Table 1) for the majority of the galaxies, both along the 
major and minor axis, at the lowest possible airmass to miti- 
gate differential refraction. We observed the Northern targets 
with APO using the Dual Imaging Spectrograph (DIS), with a 
1"5 slit, a 1200 lines mm"' grating, a central wavelength of 
4820 A, a spectral range of 4250-5400 A and a dispersion of 
finstr ~ 52 km s"', as measured from the arc lines. These ob- 
servations were performed at the galaxy major axis only. We 
also present one calcium triplet AA8498 , 8542, 8662 (CaT) mea- 
surement, based on a red setting with a central wavelength of 
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Fig. 2. — Example fits to four objects in the sample. Left: Stellar velocity dispersion fits to two B&C spectra, using the best-fit weighted sum of an A9, G6, 
G8, KO, and K3 star. We plot the B&C spectra along the major axis (middle; thick histogram), the same along the minor axis (top; thin histogram), best-fit models 
(red solid), and residuals (bottom; thin histogram), while regions excluded from the fit are demarcated with grey bars. Flux density is normalized to the continuum. 
Right: Spectra displayed as above, but we have only major axis spectra observed with DIS. Here we used an F2, G6, Kl, and K4 star 

175 U 



8540 A (8000-9 1 00 A) and comparable instrumental resolution. 

The reduction of both the DIS and B&C data proceeded in a 
similar fashion. Flat-fielding, bias-subtraction, and wavelength 
calibration were performed within irafQ In the case of the 
B&C data, flux calibration and telluric correction was also per- 
formed within iraf, using LTT 3218, LTT 3864, LTT 6248, 
Hilt 600, LTT 3218, Feige 56, and CD 32 as flux calibrator stars 
in order to match the airmass and time of observation. Flux cal- 
ibration and telluric coiTection of the DIS data used IDE rou- 
tines as described by Matheson et al. (2008). In this case we 
used Feige 34 as the flux calibrator In the Appendix we com- 
pare the [O III] luminosities between two independent data sets 
and demonstrate that our flux calibration is good to ^ 30%. 

The proper physical aperture used to calibrate the Mbh - c* 
relation remains a matter o f debate (e.g., Ferrares e & Merrita 
I2OOOI: iTremaine et al.ll2002h . and at present robust r^ measure- 
ments are only available for a small fraction of the maser sam- 
ple. Thus, we have chosen to extract spectra with apertures 
ranging from the resolution limit of the observations (^ 2") to 
the S/N limit of the spectra (^ 10"). Typically the signal-to- 
noise (S/N) ratio is optimized with apertures of ~ 4". The re- 
sulting dispersion in ct, values is included in the error budget, 
and is < 10%. 

4. STELLAR VELOCITY DISPERSIONS 

We use direct-pixel fitting (e.g., iBurbidge et al.l Il96l[) to 
measure stellar velocity dispersions. The stellar templates 
are broadened and fitted in pixel rather than Fourier space. 

^http://iraf.noao.edu/ 
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Fig. 3. — Compaiison between our measured velocity dispersions and val- 
ues from the literature (see §4. 1 for the origin of each dispersion measure- 
ment). Note that we have not included the alternate dispersion for NGC 1194 
(o-,= 184±18kms-'). 

This method is com putationally expensive relati ve to 

Fourier me thods (e.g., Tonrv & DavisI Il979l: Isimkini 119741: 
ISargent eta l. 1977; Bender 1990), but by no means prohibitive 
by modern standards. Direct pixel fit ting provides many 



benefits over Fourier techniques (e.g., iRix&Whitel 
Ivan der Marell Il994t iKelson et al.1 IJOOOt iBarth et all 
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2002 
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iBernardi et alj 120031) . For one thing, it is easier to incorpo- 
rate noise arrays directly. Of particular importance to active 
galaxies, arbitrary regions of the spectrum (such as narrow 
emission lines) may be masked in a straightforward manner 
without introducing spurious signals into the final results. 

The details of t he code used her e are described in 
iGreene & Hoi (l2006al) and lHo et alj(l2009l) . In short, each spec- 
trum is modeled as the linear combination of template stars 
T(X) that are shifted to zero velocity, broadened by a Gaussian, 
G(X), with width a and diluted by a constant [or power-law if 
necessary; C(A)]. Finally, the entire model is multiplied by a 
polynomial [P{Xy, typically third order] to account for intrinsic 
spectral differences, reddening, and errors in flux calibration: 



M(A) = P(X){[T(X) (8) G(X)] + C(X)} 



(1) 



The nonlinear Levenberg-Marquardt algorithm, as imple- 
mented by mpf it in IDL (Markwardt 2009), is used to min- 
imize x^ in pixel space between the galaxy and the model. We 
have a library of template stars observed with the B&C spec- 
trograph. For the fitting we use the following stars: HD 36634 
(G6 III), HD 129505 (G8 IIII), HD 126571 (KO III), HIP 70788 
(K3 III), and HD 131903 (A9 V). For the APO observations, 
our spectral library is much more limited, and we have only 
one reliable velocity template star. Thus, we utilize a large li- 
brary of high-S/N, high-spectral resolut ion (^ 26 km s"') stel- 
lar templates from IValdes et al.l (|2004|) . Specifically, we use 
the following stars: HD 26574 (F2 II-III), HD 107950 (G6 III), 
HD 18322 (Kl III), and HD 131507 (K4 III). In order to trans- 
late the measured dispersion to the true dispersion, we first 
measure the instrumental dispersion of the DIS spectra using 
arc lines (52 km s"'), which yields a quadrature difference of 
(52^-26^)°^ = 45 km s"' between the DIS and Valdes spectra. 
We also measure the dispersion of the DIS velocity template us- 
ing a Valdes star of the same spectral type and find 45 km s"'. 
The consistency of the two methods is reassuring. Thus, the fi- 
nal dispersion is simply cr» = ((7,^^35 -45^)*''^. Note that because 
we use stellar templates observed with an identical instrumental 
set-up, we need apply no corrections to the B&C spectra. 

Given our spectral coverage, the most robust spectral fea- 
tures with high enough equivalent width (EW) for dispersion 
measurements are either the Mg lb triplet or the G-band. We 
avoid the Ca Hh-K features because they are highly sensitive to 
spectral type, occur in a spectral region with a strong continuum 
discontinuity, and may be contaminate d by interstellar absorp- 
tion (see review in lGreene & Hol2006ah . The Mg lb features are 
in many cases contaminated by emission from [Fe VII]A5158 A 
and [N I] AA5 198, 5200 A. Furthermore, differences in the a/Fe 
ratios between the template stars (from the solar neighborhood) 
and external galaxies will le ad to system atic errors in a^, mea- 
surements from this region dBarth et al] |2002). For these rea- 
sons, we focus on the G-band region for the velocity disper- 
sion measurements. We have multiple slit positions for most of 
the Dupont spectra. In Table 2 we present measurements for 
each slit position with the 4" aperture, but we also present the 
weighted mean of the two, and use that value in all subsequent 
work. Example fits are shown in Figure 2 while the dispersions 
are presented in Table 2. 

4.1. Literature Comparison 

There are literature measurements available for a limited 
number of targets in our sample. Direct comparison can be 



tricky; the dispersions may be highly ap erture dependent (e.g., 
iPizzeUa et aTlbOO^lBarbosa et al.ll2006l) . Nevertheless, in gen- 
eral we find decent agreement with the literature values, as 
summarized in Figure 3. Two tai-gets (NGC 2273, NGC 4388) 
appear in the recent atlas of stellar velocity dispersions mea- 
sur ed for the Palo mar spectroscopic survey of nearby galax- 
iesj Ho et al. 19951) . taken with a comparable 2" x 4" aperture 
JHo et al.l l2009). The two sets of observations agree within the 
uncertainties. For NGC 2273 we find 146 ±17 km s"' and Ho et 
al. present 149 ±9 km s"', while for NGC 4388 the values are 
104 ± 10 km s"' along the major axis and 92 ± 10 km s"' along 
the minor axis. In the case of NGC 4388, note that our observa- 
tion is only marginally spectrally resolved. However, our SDSS 
measurement, presented in the Appendix, is 107 ± 1 1, giving us 
some confidence in this result. Oth er literature co mparisons for 
these two targets are presented in Ho et al.l (|2009|) . 

Two other galaxies (N GC 3393 and IC 2560) were observed 
bv lCid Fernandes et al.l (f2004) with a l."5 x 0."82 slit. Again, 
the measurements agree within the quoted errors. In the case 
of IC 2560, we find 134± 12 km s~' along the major axis com- 
pared to their 144 ±20 km s"', while for NGC 3393 we find 
142 ± 16 km s"' along the major axis and 156 ± 16 km s"' along 
the minor axis, as compared to 157 ± 20 km s"'. NGC 3393 is 
also presented bv lTerlevich et al.l d 19901) measured within a 2" 
extraction. They find 1 84 ± 1 8 km s"' , which is marginally con- 
sistent with our minor-axis measurement. 



4.2. Emission Line Spectra and Bolometric Luminosities 

One convenient by-product of our ct, measurement technique 
is a pure emission-line spectrum that has been corrected for un- 
derlying absorption lines (e.g., H/3). Since the narrow emission 
lines (particularly [O III]) represent one of the only methods 
for estimating bolometric luminosities for these obscured ob- 
jects, we descri be the measurements here . The fitting routine s 
are described in lGreene &"Ho! (2005) and lGreene et al.1 (|2009|) . 
Briefly, we fit the [O III] A4959,5007 lines with sums of up 
to four Gaussians. The two lines are restricted to have identi- 
cal shapes and the velocity difference and relative strengths are 
fixed to laboratory values (3.1 in the latter case). The H/3 line 
is fitted independently with the sum of two Gaussians (which 
is always completely sufficient given the typical signal-to-noise 
ratio, S/N, in this line). We also calculate the full-width at half 
maximum (FWHM) of the summation of all components in the 
[O III] model. The velocity offsets between individual compo- 
nents are small compared to the total line width, so the FWHM 
is well-defined. To bracket systematic uncertainties, we mea- 
sure the line width in the spectra from the 4"and 10" apertures 
and the error bar represents the difference in those two widths. 

For the majority of the B&C spectra, as well as the DIS ob- 
servation of NGC 6323, we also model the HaH-[N II] com- 
plex. Ideally we would use measurements of the [ S II] Une 
width as a model to debl end Ha from [N II]. (e.g., iHo et al.l 
1 19971; iGreene & Ho! '2005). However, in this case we simply 
describe each narrow line with the same combination of Gaus- 
sians, fix flie relative strengths of the [N II] AA6548,6583 A 
lines to 1 :2.96, and fix the relative wavelengths of all three lines 
to their laboratory values. In the case of NGC 4388, we fur- 
ther include a b road com ponent fixed to have the same FWHM 
as reported bv .Ho et al.l (11997 ). but we are not sensitive to the 
presence of weak broad lines in the other objects. We derive an 
intrinsic reddening correction for [O III] using the H/3 and Ha 
fluxes when available, but do not apply it to the fluxes. All the 
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Table 2. Galaxy Properties 
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logo-. 


log tT„,„i 


log CTmin 


r 


Ar 


A, 


B/T 


g-r 


T,. 


log Mbuigc 


log Mbh 


,/lGyr 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


NGC1194 


2.17±0.07 






14.0±0.3 


0.21 


0.29 


0.5±0.2 


0.7±0.3 


2.6±1.6 


10.3±0.4 


7.82± 0.05 


< 0.05 


NGC2273 


2.16±0.05 




















6.88± 0.05 




UGC3789 


2.03±0.05 




















7.05± 0.05 




NGC2960 


2.22±0.04 


2.17±0.04 


2.26±0.04 


14.3±0.1 


0.12 


0.17 


0.4±0.1 


0.5±0.1 


2.0±0.1 


10.2±0.1 


7.05± 0.05 


0.19±0.13 


IC2560 


2.15±0.03 


2.13±0.03 


2.18±0.03 


















0.20±0.13 


NGC3393 


2.17±0.03 


2.15±0.03 


2.19±0.03 
















7.49±0.12 


0.()2±0.01 


NGC4388 


2.03±0.03 


2.02±0.03 


2.05±0.03 


12.0±0.4 


0.09 


0.13 


0.5±0.2 


0.4±0.4 


1.3±0.1 


9.79±0.3 


6.93± 0.05 


<0.05 


NGC6264 


2.20±0.04 


2.22±0.04 


2.06±0.04 


15.4±0.2 


0.18 


0.25 


0.5±0.1 


0.4±0.2 


1.4±0.1 


10.2±0.1 


7.45± 0.05 


O.OOliO.OOl 


NGC6323 


2.20±0.07 






16.0±0.1 


0.05 


0.07 


0.2±0.1 


0.7±0.1 


2.8±0.1 


lO.OiO.l 


6.96± 0.05 





Note. — Col. (1); Galaxy. Col. (2): ct* within ~ 4" (km s"'); weighted average of two slit positions when available. Col. (3); a^ within ^ 4" (km s"^) measured 
along the major axis. Col. (4): cr* within ^ 4" (km s"^) measured along the minor axis, when available. Col. (5); r-band bulge magnitude (mag), as measured from 
SDSS images using GALFIT. Col. (6): Galactic reddening (mag) in the r-band. Col. (7); Galactic reddening (mag) in the g-band. Col. (8); Bulge-to-total ratio in the 
r-band. Col. (9): g-r color (mag) measured from SDSS images, fixing the g-band model to that derived from r. Col. (10); Mass-to-light ratio (T©) of the bulge in 
the r-band inferred from the g-r color using the relations presented in Bell et al. 2003. Note that the values are considerably lower than in elliptical galaxies. Col. (1 1): 
Bulge mass (Mq) as estimated from the r-band luminosity (Col. 5) and the Tr (Col. 10). Col. (12): Black hole mass (Mq) as measured using megamaser spots that are 
in Keplerian rotation fractions of a pc from the black hole (Kuo et al. in prep.). Col. (13): Fraction of light in a 1 Gyr population as estimated from fits to our spectra using 
Bruzual & Chariot (2003) stellar population models with an exponennal decay dme of r = 1 Gyr. The spectra are fit with two components, one with an age of 1 Gyr. the 
other with an old population of 10 Gyr. In all cases the 10 Gyr population dominates the light in the aperture. 



line measurements are summarized in Table 3. 

We are primarily interested in the [O III] luminosities as an 
isotropic indicator of the AGN luminosity. Estimating robust 
bolometric luminosities in these objects is complicated since 
megamaser sy stems are fou nd preferentially in obscured active 
galaxies (e.g. jBraatz et al.l il997). Among the syst ems with cir- 
cumnuclear disks, > 75% are Compton thick (e .g., IZhang et aH 
l2006HGreenhill et al.ll2008tlZhang et al.ll2010l) . Thus, there are 
only a few available luminosity diagnostics that are thought 
to be isotropic and immune from significant absorption. One 
is the strong and ubiquitous [O III] line that has been used 
both t o find AGNs and to determine their intrinsic lu minosities 
(e.g., iKauffmann et al.ll2003t IZakamska e t al.ll20 03|) although 
concerns about redde ning remain (e.g.. iMulchaey etaDll994t 
iMelendez et al.l 120081) . The shape and equivalent width (EW) 
of the Fe Ka line at 6.7 keV are also sensitive to the intrin- 
sic luminosity and in clination of the obscuring material (e.g., 
iLevenson et al.ll2002h . In a few cases, detailed modeling of the 
Fe Ka lines have been performed, which yie lds an indepen- 
den t estimate of the bolometric luminositie s (ILevenson et al.l 
12006) . Finally, iDiamond-Stanic et al.1 (l2009h argue for the use 
of the [O IV] A25.89 /xm emission Une as a higher-fidelity tracer 
of AGN bol ometric luminosity (see also IMelendez et al.ll2008t 
iRigby et al.ir 2009). 

Using the [O III] bolometric correction presented in lLiu et al] 
(I2009.) . we derive bolometric luminosities for the entire sample 
using our measurements of L[om] (Table 3). Since this bolo- 
metric correction was derived without correction for internal 
reddening, we do not deredden either, but tabulate the proper 
correction for reference and comparison with the other indi- 
cators. The estimated Eddington ratios range from Lboi/iEdd= 
10~^-0.16. Including internal reddening would change the val- 
ues by 0.6 dex at most. There are only two objects with esti- 
mated bolometric luminosities from the Fe Ka line (excluding 
IC 2560 for the moment), also shown in Table 3. There are 
also three cases with [O IV] measurements in the literature; we 
have converted to bolometric luminosity using the calibration of 
J. Rigby et al. in preparation. As expected, Lboi([0 IV]) tends 
to be higher than Lboi([0 III]), but the number of objects is too 
small to draw definite conclusions about the average magnitude 



of this effect, or how directly it is related to our internal redden- 
ing estimates. 

4.3. Stellar Populations 

The presence or absence of young stellar populations in the 
bulge regions of these galaxies provide important clues as to 
their growth and evolution (§5.3). The spectra provide useful 
constraints on stellar ages, although the modeling is compli- 
cated by the narrow emission lines that litter the spectrum. We 
do not perform full stellar population synt hesis modeling. In- 
stead, we take a few representative sets of JBruzual & CharlotI 
(2003) models and perform a least-squares fit to the galaxy con- 
tinua, fixing the stellar velocity dispersion to that derived a bove. 
We adopt the Pa dova isochrones (Marigo & Grrardil2007l) and a 
IChabried (l2003l) initial mass function. We take both single-age 
stellar population models and exponentially decaying models 
with a decay time r = 1 Gyr. For simplicity, we use two models 
with ages of one and 10 Gyr respectively, although we experi- 
ment with models having ages of 0.5 and 5 Gyr and also one- 
half of solar metallicity. We adopt a singl e reddening sc r een fo r 
each galaxy, using the reddening law of ICalzetti et al.l (|2000|) . 
See example fits in Figure 4. The fraction of light in young 
stars is tabulated in Table 2. 

Of the six galaxies for which our spectra include the 4000A 
break, we see that only two show evidence for a significant con- 
tribution from intermediate-age populations (NGC 2960 and IC 
2560). NGC 3393 and NGC 6264 may have a very minor young 
component, but presumably higher spatial resolution would be 
needed to isolate it. We can compa re these results with t h e mor e 
detailed stellar population fits of ICid Fernandes et al.l ( |2004|) . 
whose sample also includes IC 2560, NGC 4388, and NGC 
3393. In NGC 3393 we both agree that the stellar populations 
are old; Cid Fernandes et al. find that 4% of the luminosity in 
this source come from stars with an age < 25 Myr. They find 
a substantial 33% contribution from young stars in IC 2560, 
consistent with our findings. The only ambiguous case is NGC 
4388. They find that 13% of the luminosity comes from this 
young component, while we do not find compelling evidence 
for a young component. In fact, NGC 4388 is quite difficult to 
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Galaxy /h^ fuc. /loilli 

(1) (2) (3) (4) 



'[OIII] 

(5) 



(6) 



log i|OI 
(7) 



(8) 



log A/bh 
(9) 



I ^bol,IOIII]^^Etld log /-boLx^^Edd 'og Abul,IOIV]^^Edd 
(10) (II) (12) 



NGCII94 


0.09 




30 




430± 43 


39.9 




7.82 


-2.17 


NGC2273 


0.14 


0.75 


380 


5 


190± 19 


40.5 


40.1 


6.88 


-1.64 


UGC3789 


0.08 


0.47 


700 


3 


200± 23 


41.3 




7.05 


-0.82 


NGC2960 


0.09 


0.31 


50 


1 


700± 70 


40.5 




7.05 


-1.65 


IC2560 


0.08 


0.26 


450 


2 


360± 36 


41.0 


41.1 






NGC3393 


0.10 


0.24 


630 


1 


460± 46 


41.3 




7.49 


-0.80 


NGC4388 


0.08 


0.29 


530 


2 


380± 38 


40.4 


41.6 


6.93 


-1.76 


NGC6264 


0.10 




80 




550± 55 


41.3 




7.45 


-0.88 


NGC6323 


0.28 




20 




740±255 


40.4 




6.96 


-1.75 



-2.85 



-2.06 



0.42 



Note. — Col. (1): Galaxy name. Col. (2): Rux in H/3 normalized to flux in [O ill] (Col. 4). Corrected for Galactic reddening only (see Table 2). Col. (3): Ha 
flux normalized to [O ill] flux (Col. 4). Corrected for Galactic reddening as above. Col. (4): Flux in [O in] A5007 A line (10"'^ erg s"' cm"-). Corrected for Galactic 
reddening as above. Col. (5): Internal reddening coiTcction to ]0 in] as derived using the Balmer decrement (Cols. 2, 3). We have not apphed this correction to 
any emission-line measurements. Col. (6): FWHM in the ]0 ill] line (km s"^) as measured within a 4" aperture to match the cr* measurement, with the instrumental 
resolution subtracted in quadrature. Col. (7); Luminosity (erg s"') in the ]0 in] line. Col. (8): Luminosity in the [O iv] A25.89 /im line from Diamond-Stanic et al. 
2009. Col. (9): Black hole mass (Mq) as measured using megamaser spots from Kuo et al. 2010 and Ishihara et al. 2001. Repeated here for convenience. Col. (10): 
Inferred Eddington ratio, Lboi/LEdd with the bolometric luminosity (erg s"^) based on Ljoiii] and the bolometric coiTcction from Liu et al. 2009. Col. (1 1): Inferred 
Eddington ratio, Lbol/^Edd with the bolometric luminosity (erg s"^) based on modeling of the Fe Ka line from Levenson et al. (2006). Col. (12): Inferred Eddington 
ratio, Lbol/^Edd with the bolometric luminosity (erg s"^) based on the measured Ljoivj (Col. 8) and the bolometric correction of Rigby et al. in preparation. 



model, both because it is extincted and because the emission 
lines have very high equivalent widths. 

Those galaxies (NGC 2273, UGC 3789, and NGC 6323) that 
we observed with APO do not have the blue coverage needed 
to constrain the stellar populations. In the case of NGC 2273, 
young stars asso ciated with an inner ring are well-documented 
dGu et alj|2003h . Unfortunately, we do not have much infor- 
mation for the other two. The bulge color of NGC 6323, 
r = 0.7 (see below ), is typical for an early-type spiral galaxy 
Fukupta et alJI 19951) . suggesting that either there is not much 
star formation and/or that there is a considerable quantity of 
dust. Overall we find evidence for ongoing or recent star for- 
mation in roughly half of the sample. 

4.4. Uncertainties: Aperture Effects, Rotation, and Bars 

As discussed briefly in §2, the velocity dispersion pro- 
files of late-type spirals can be complicated by kinematically 
cold, disk-like components as well as bars or other non- 
axisymme tries. As a result, it is impossible to pick a single, 
well-motivated aperture size, and the long-slit data presented 
here are not deep enough to derive robust rotation curves. We 
can say that the measured dispersions do not change by more 
than 10% over the range of radii we probe, and this spread in 
values is included in the error budget. Robust r^ measurements 
are needed to define the optimal radius of extraction for each 
galaxy. Along similar lines, it is possible that the s lit orienta- 
tion r e lative to that of the bar might impact a^, (e.g.. lKormendvl 
119831; lEmsellem et alj|2001h . However, for the objects with 
both minor and major axis spectra, we see at most a 40% differ- 
ence in the two measurements. Ideally, we would like to obtain 
integral-field spectroscopy for the objects to settle this ambigu- 
ity. 

An additional concern is that differing stellar populations and 
metallicities lead to bias due to template mismatch. While our 
code includes a mix of spectral types, subtle s ystematic s from 
region to region remain (see an example in iHo et al.ll2009h . 
When we calculate errors in cr*, we include the differences be- 
tween the nominal dispersions and those measured from the 
5100-5400 A region [excluding the Mg lb lines themselves 



since variations in ct/Fe tend to b ias the dispersion measure- 
ments high (e.g.. lBarthetal.ll2002l) 1. 



5. PHOTOIVIETRIC DECOIVIPOSITION 

Often, spiral galaxies contain bars, ovals, rings, dust lanes, 
nuclear spirals, and nuclear st ar cluster s in their centers (e.g., 
ICarollo et al.ll997HBoker eFal . 2002; Martini et al.l2003b . Fur- 
thermore, in galaxies with young stellar populations, the mass- 
to-light ratio (T) may vary significantly both within a given 
galaxy and from object to object. To properly decompose such 
light profiles requires both deep imaging on large scales to 
quantify the disk component and high resolution nuclear imag- 
ing, ideally in a red band, to reveal the complex nuclear mor- 
phology. 

We first tried to perform detailed profile decomposi- 
tions using Two-Mic ron All-Sky Survey (2MASS) images 
JSkrutskie et al.ll2006h . As shown in Figure 5, these images are 
too shallow to detect large-scale disks and rings. Furthermore, 
the resolution (2-3") is insufficient to uncover nuclear features. 
A combination of TiT-band observations with HST and from the 
ground is needed to fully model the entire galaxy. As an initial 
exploration we have settled on a compromise. There are five 
galaxies that fall within the Sloan Digital Sky Survey footprint. 
We have performed two-dimensional image decomposition on 
the gr images of these five galaxies, which allows us to measure 
a magnitude and a color for their bu lge components. We then 
use the relations of iBell et al.l (l2003h to convert the g- r color 
into Tr and thereby stellar mass. 



5.1. GALFIT 

We perform image decomp osition in two dimensions using 
the program GALFIT dPeng et al. 2002). Two-dimensional fit- 
ting allows flexibility in the ellipticities and position angles 
of different components that can hel p break degeneracies in 



fitting multi-coniponen t galaxies (e.g., lAndredakis et alj 19951; 
IWad adekareTan[T999h 
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as reviewed in detail in iPeng et akl 
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Fig. 4. — Three examples of fits to the entire spectral range using high- 
resolution Bruzual & Chariot 1 2003) models (see §4.3). We plot the oiiginal 
spectra (thick black histogram), the best-fit model with cr, fixed {red solid) and 
the residuals (thin black histogram). Flux density is arbitrarily normalized and 
masked regions are indicated with grey bars. 

GALFIT and related programs create model galaxies as the 
sum of ellipsoids with parametrized surface brightness profiles 
that are convolved with the point-spread-function (PSF) of the 
image and then compared with the data. Best-fit parameters are 
those that minimize the differences between data a nd model. In 
general, we model the galaxies with lSersig ( 1 19681) functions: 



sc loa 




1 5 10 

Fig. 5. — SDSS r-band images (left) and radial profiles (right) for the 
galaxies NGC 1194 (top) and NGC 6264 (bottom). A scale of 10" is shown. 
The radial profiles show the S'-band (open circles) and SDSS r-band (filled 
circles) isophotes, in order to demonstrate that the optical data are both deeper 
and have higher resolution. 

An exponenti al disk has n = 1 whil e n = 4 corresponds to 
the well-known Ide Vaucouleu rs (1948) profile. We adopt the 
PSF models provided by the SDSSQ, and note that our analysis 
is not very sensitive to the detailed PSF shape. Unlike typical 
two-dimensional modeling, GALFIT now incorporates coordi- 
nate rotations and Fourier modes (departures from axisymmet- 
ric modeling) that are particularly useful in handling dust-lanes, 
bars, and spiral structure (Peng et al. 20101) . 

As described above, in addition to the structural measure- 
ments we also need color information to estimate the stellar 
masses of the galaxies. We first fit the r-band images to ex- 
tract structural information and then apply the same model to 
the g-band images, allowing only the total flux to vary for each 
component. In Figure 6 we show the fits to the r-band images. 
We caution that the measurements are preliminary in the ab- 
sence of high-resolution imaging. Therefore, we present only 
the "bulge" luminosities and colors, and reserve more detailed 
discussion of the galaxy structure (e.g., Sersic index) to future 
work. For the moment, we simply note that the bulge colors 
provide an additional crude handle on the stellar populations in 
these galaxies (Table 2). In the case of NGC 1194 and NGC 
6323, the g-r color is ty pical of Sa galaxies dColeman et alJ 
ll980l:lFukugita et alJ[T995 ). NGC 2960 has the color of an Sbc 
galaxy, and NGC 4388 and NGC 6264 are the bluest, with col- 
ors typical of Scd galaxies. 



/(r) = 4 exp { -b„ 



\/n 



(2) 



where r^ is the effective (half-light) radius, I^ is the intensity at 
r^,, n is the Sersic index, and b„ is chosen such that 



r- 

I{r)2'Krdr = 2 / I{r)2iirdr. 
10 Jo 

■'http://www.sdss.org/dr7/products/images/read_psf.html 



(3) 



5.2. Uncertainties in Decomposition 

Clearly, errors in the PSF shape and sky level add scatter to 
the photometric decompositions. We expect that these uncer- 
tainties are < 0.1 mag. Therefore, we are dominated by system- 
atic errors. More specifically, because of prominent dust lanes 
in many galaxies, combined with the limited spatial resolution 
of the SDSS, we do not have a strong constraint on the Sersic 
index of the bulge, which in turn contributes large errors to the 
effective radii as well as in the ratio of bulge-to-total luminosity 
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Fig. 6. — Photometric fits to all objects with SDSS images. We show the r-band image (left) with a scale-bar of 20", the best-fit GALFIT model (middle), and the 
residuals (right), as well as radial profiles for the data (black circles), the total model (solid black line), any bulge (red solid line), disk (blue solid line), bar (green 
solid line), or outer disk (magenta solid line). Residuals are shown beneath in solid squares. Note that dust-lanes and spkal arms can infi'oduce non-monotonic 
behavior in the radial profiles. 
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Fig. 7. — Three im ages of NGC 2273 at progressively smaller scales. The lai'ge-scale i?-band images are from WIYN and were kindly provided by P. Erwin 
jErwin & Sparkel2003l ; left, middle), while the nuclear spiral is a NICMOS/F160W image. Finally, we show the radial profiles for the NICMOS (blue open circles; 
top), WFPC2/F606W {red filled circles), and i?— band (crosses) images. 



(B/T). 

Currently, our error bars represent the spread in derived lu- 
minosities resulting from a range of different models, which 
give some indication of the uncertainties due to model assump- 
tions (e.g., Peng et al. 2010). After running these models, we 
report the one with the lowest value of reduced y^ and assign 
parameter uncertainties based on the spread across the models. 
Let us take as an example the galaxy NGC 1194. In this case 
we have fit a model with the bulge Sersic index fixed to n = 4 
and the disk fixed to « = 1 . As the next level of complexity we 
have allowed each component to vary freely, yielding n = 3 .9 for 
the bulge and n = 2.2 for the "disk". Finally, we have run two 
models with different parameterizations of the dust lane that 
have n w 2.7,« w 1.5 for the bulge and disk respectively. In all 
cases the fits would be deemed acceptable, but using a range of 
models allows us to quantify the dependence of the bulge pa- 
rameters on our model assumptions. In this case, the resulting 
bulge magnitude varies by ^ 0.5 mag between the model with 
n = 4 and those including dust. These systematic uncertainties 
exceed by far those resulting from sky and PSF uncertainties, 
and are the ones we report here. 

5.3. Bulge Classification: Classical or Pseudobulge 

In this subsection we discuss the morphological properties 
of the maser sample. In particular, we address the bulge mor- 
phology. First consider NGC 1194 (Fig. 5a). This galaxy has a 
massive bulge component and is classified as an SO. Both in the 
radial profile and the image, it is clear that the bulge in NGC 
1194 has smooth isophotes and a large ratio of bulge-to-total 
(B/T) luminosities. Much more typical of the megamaser disk 
galaxies is NGC 6264, which has a strong bar (Fig. 5b). Other 
barred galaxies include NGC 2273, NGC 3393, NGC 4388, and 
IC 2560. Roughly speaking, the bar fraction in the maser sam- 
ple is consist ent with the ^ 60% seen in the spiral galaxy pop- 
ulatio n (e.g.. lEskridge et all 120001 : iMenendez-Delmestre et alj 
120071) . 

Finally, we examine NGC 2273 (Fig. 7), also a barred spi- 
ral galaxy. In addition, NGC 2273 contains two outer rings 
and a centrally concentrated bulge component dErwin & Sparkd 
l2003h . Zooming in with HST/WFPC2 and NICMOS reveals 
that the galaxy center consists of an inner ring surrounding 
a nu clear spiral (e.g.. iMulchaev et al.l 119971 : lErwin & Sparkd 
l2003b . This spiral component is very different in shape, kine- 
matics and stellar populations from a classical bulge. The latter 
is understood to be a dynamically hot stellar component with 
an old stellar population, effectively an elliptical gala xy with 
a surrounding disk (see discussion in .Kormendv & KennicuttI 



|2004|) . In contrast, the center of NGC 2273 is a disk with ongo- 
ing star formation fueled by a reservoir of molecular gas (e.g., 
iPetitpas & Wilson! I2OO2I) that is likely fed by the larger-scale 
bar Furthermore, there is a distinct drop in a^, within the inner 
5" where the kinematics are dominated by rotation, presumably 
because the young bright stars trace the motion of the gas in 
which they formed (iBarbosa et al.ll2006t iFalcon-Barroso et al] 
2006). 

It has been clear for a long time that the "bulges" 
of spiral galaxies are often disk-like. The first clues 
were kinematic, with some bulge-like components display- 
ing very high ratios of rot a tional velocity to velocit y dis- 



. „ |ty 

persion (e.g., BinnevI 119781 : iKormendv & IllinpworthI 119821 : 
lKormendvlll993 '). HST imagi ng reveals disk-l i ke and flattened 
centra l isophotes (e.g.. iFathi & Peletieij 120031 iF isher & Droryl 
2008; 



Andredakis & SandersI 
l2003lv 



19941: 



Courtea u et al.l 
nuclear spirals 

(e.g 



19961: 



iMac Arthur et al.l 120031) . bars, ovals, nuclear spi rals and 
rings in many bulge-lik e components (e.g., Car oUo et al.l 
[1997; Martini e t al.l 12003 ) as well as ongoing star formation 
(iFisher et al...2009l ). To distinguish these systems from "classi- 
cal" bulges, they are referred to as "pseudobulges". It is thought 
that secular, internal processes (e.g., bars or spiral arms) are 
responsible for building pseudobulge mass, in contrast to the 
rapid merging processes, i.e. violent relaxation, thought to 
build elliptical gala xies and classical bulges . The literature is 
reviewed in detail in lKormendv & KennicuttI ( |2004|) . Through- 
out, we will use the term "bulge" to refer to any luminous, cen- 
trally concentrated galaxy center, where classical bulges are of 
the elliptical type and pseudobulges are essentially disks. 

Different author s use different (overlapp ing) criteria to iden- 
tify pseudobulges. IFisher & Drorvl (|2008|) start with a sample 
of morphologically selected pseudobulges, meaning that a nu- 
clea r bar, spiral, or ring is detected within the bulge region (see 
also iKormendv & Kennicutal2004l) . They show that the distri- 
bution of Sersic indices in local bulges is bimodal, and sug- 
gest that pseudobulges are galaxies with « < 2 (for a defini- 
tion of the Sersic index, see §5.1). These authors also show 
that the average B/T of pseudobulges (0.16) is lower than that 
of classical bulges (0.4) with a large s pread. iGa dotti (2002), 
on the other hand, advocates use of the Kormendvl (11977) rela- 
tion as a discriminator, since pseudobulges tend to ha ve lower 
central surface brightnesses at a fixed radius (see also ICarollol 
1 19991: IFisher & Drorvl2008h . Finally, while classical bulges are 
typified by old stellar populati ons, pseudobulges tend to have 
ongoi ng star f ormation (e.g., Kormendv & Kennicujd 120041: 
Drorv & Fisheii!2007: Fisher et al. 2009: Gadotti 20091k Since 
deriving robust velocity measurements is beyond the scope of 
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this paper, and in the absence of more robust structural informa- 
tion, we rely on morphology and stellar population properties 
at the present time. 

The nearest, well-studied targets in our sample (NGC 4388, 
and NGC 2273) probably contain pseudobulges. In the case of 
NGC 2273 this classification is based on both the young stel- 
lar populations and the rings and nuclear disk. NGC 4388 is 
less certain, but there is clear evidence for recent star formation 
and dust. We suspect that NGC 6264 contains a pseudobulge, 
given its morphological similarities with NGC 2273 (namely 
the outer ring and inner bar) and the evidence for young stars. 
The same goes for NGC 3393 and IC 2560, which each contain 
an outer ring, a bar, and an inner ring. On the other hand, NGC 
1194, with both evolved stellar populations and a large bulge, 
probably contains a classical bulge. NGC 2960 has some of the 
clearest evidence for ongoing star formation, and so we tenta- 
tively put it into the pseudobulge category. Finally, we remain 
agnostic about NGC 6323, which is one of the most distant tar- 
gets. Thus, of the nine targets we consider, at least seven likely 
contain pseudobulges. 

6. SCALING BETWEEN Mbh AND 0-* 

In Figure 8 we present the location of the megamaser galax- 
ies in the Mbh - f* plane. The maser galaxies do not fol- 
low the extrapolation of the Mbh - cr* relation defined by the 
elliptical galaxies. Instead, they scatter towards smaller BH 
masses at a given velocity dispersion. Quantitatively, tak- 
ing A MBH=log(MBH) - log[M(cr*)], where log[M(cr*)] is 
the expected Mbh given a^, we find (Z\Mbh) = 0.24 ±0.10 
dex. There are many hints in the literature that the Mbh- 
(T* relation does not extend to low-m ass and la t e-type galax- 
ies in a straightforward manner (e.g., lHull2008t iGreene et al.l 
[2008; Gadotti & Kauffmann 2001). However, the precision BH 
masses afforded by the maser galaxies make a much stronger 
case. The Mbh - c* relation is not universal. Neither the shape 
nor the scatter of the elliptical galaxy Mbh - cr* relation pro- 
vides a good description of the maser galaxies in this plane. 

We now add the maser galaxies to the larger sample of lo- 
cal galaxies with dynamical BH masses to show that indeed a 
single, low-scatter power-law does not provide an adequate de- 
scription of all galaxies in the Mbh - cr* plane. For convenience 
and to facilitate comparison with previous work, we assume a 
power law for all fits, although that form may not provide the 
best description of the sample as a whole. 

6.1. Subsamples 

We fit not only the full sample, but also divide the galaxies 
by morphological type. For ease of comparison, we adopt the 
divisions of Giiltekin et al. (2009). In essence, the galaxies are 
divided into two primary morphological bins: elliptical and spi- 
ral galaxies. For completeness we perform fits assuming that SO 
galaxies belong to each subgroup. We adopt nomenclature from 
Giiltekin et al.. "Elliptical" refers only to the elliptical galaxies, 
while "early-type" refers to elliptical and SO galaxies. "Late- 
type" galaxies are spiral galaxies excluding SOs while "nonel- 
liptical" includes all spiral galaxies (SO and later). All fits are 
shown in Table 4 and Figure 9. 

In addition to these two divisions, we also consider pseudob- 
ulges by themselves. Early work suggested that pseudobulges 
(see §5.3) obey the Mbh - cr* relation of class ical bulges based 
on a very small sample of galaxies ( Kormendv & Gebhardj 
l200lh . Further study, based on larger samples, has pointed to 



differences both in the slope and scatter of BH-bulge relat ions 
in ps e udobulges JHu 2008; Greene et al. 2008; Graham &Lil 
120091; iGultekin et all l2009t iGadotti & Kauffmann 2009), but 
the results are not conclusive, predominantly because the avail- 
able number of dynamical BH masses in pseudobulges is so 
small. The maser sample mitigates this problem. As a fidu- 
cial sample, we adopt the pseudobulges tabulated by Giiltekin 
et al.. However, both NGC 3245 and NGC 4342 were classified 
as classical bulges by Fisher & Drorv (2008) and Kormendv! 
(2001) respectively. Thus, we have performed a second fit 
(Pseudo2) excluding these two galaxies. While the slopes differ 
between the two fits (although not significantly), the intrinsic 
scatter does not change. 

6.2. The Fitting 

To facilitate direct comparison, we adopt the fitting proce- 
dure of Giiltekin et al.. As is commonly done, we fit the rela- 
tion: log (Mbh/M0)= ;9 + alog(cr*/200km s"'). We use a gen- 
eralized maximum likelihood analysis that is designed to in- 
corporate intrinsic scatter in a natural manner For simplicity, 
we assume a Gaussian distribution in both the measurement er- 
rors and the intrinsic scatter (directly investigated by Giiltekin 
et al.). For a set of observed points (M,, cr*,), we maximize the 
total likelihood: 

/: = []z,(M,,a,,,). (4) 

In the presence of measurement errors, if the likelihood of mea- 
suring a mass Mobs for a true mass M is Qi{Mobs\M)dMob^, then 
for a given observation the likelihood is: 



li= / Qi{Mi\M)P{M\a,MM 



(5) 



We assume both Q and P to have a log-normal form. Upper 
limits are included, where the mass is taken to be below the up- 
per limit with a Gaussian probability given by the measurement 
uncertainty (see details in Appendix A of Giiltekin et al.). Their 
inclusion does not change the fits significantly. 

Since errors in the independent variable are not naturally in- 
cluded, we run 2000 Monte Carlo simulations for each fit in 
which the dispersion values are drawn from a Gaussian distri- 
bution around the measured value, assuming symmetrical errors 
in the logarithm. The reported error is a quadrature sum of the 
error from the fit and the width of the final distribution of fit- 
ted values from the Monte Carlo runs. As a s anity check, we 
also em ploy a x^ fitting technique adapted from lTremaine et aU 
(|2002). In this formalism, intrinsic scatter is included as an 
additional term such that the final reduced %-= 1. While we 
only present results from the maximum likelihood treatment, 
we note that the fits are consistent between the two methods in 
all cases. 

The formal fitting quantifies the visual impression of Figure 
8. At low BH mass (MbhJS 10^ -^o)^ galaxies begin to deviate 
significantly from the Mbh - cr, relation of elliptical galaxies. 
The zeropoint is lower by 0.4 ± 0.2 dex and the scatter is larger 
by 0.2 ± 0. 1 dex. It is worth noting that, in agreement with Giil- 
tekin et al., the fits to the pseudobulge galaxies are in no way 
distinguishable from fits to either the late-type or nonelliptical 
subsamples. Our data demonstrate that the power-law relation- 
ship between Mbh and cr, is not universal. Specifically, we find 
a population of BHs that apparently have not fully grown to the 
present day Mbh - cr* relation of elUptical galaxies. 

6.3. Potential Caveats 
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Fig. 8. — The relation between BH mass and bulge velocity dispersion for the maser galaxies presented here (open circles) and those from the literature (grey 
stars). IC 2560 is in dicated with a cross and the BH mass eiTor bar is heuristic only. For reference, we show the Mbh ~c* relation of elliptical galaxies from 
IGultekin et alj 120091 red dashed line). The maser galaxies trace a population of low-mass systems whose BHs lie below the Mgn ~ f* relation defined by elliptical 
galaxies. The largest outlier galaxies are (from highest to lowest Mbh) NGC 2960, NGC 6323, and NGC 2273. 



Our measurement techniques may not yield a^^ values that 
are directly comparable with those in elliptical galaxies. There 
is firstly the question of whether we are contaminated by disk 
light. From the SDSS imaging, it is quite clear that the bulge 
completely dominates the light within 5". As an additional san- 
ity check, in Figure lOa we look for correlations between the 
galaxy axial ratio (as proxy for the inclination) and deviations 
from the Mbh-ct* relation, A Mbh. There is no compelling ev- 
idence for a correlation between the two (Kendall's tau value 
of T = -0.2 with a probability P = 0.7 that no correlation is 
present), suggesting that our cr* measurements are not strongly 
biased by rotation in the central regions. 

We still worry about aperture biases in our ct* measurements, 
because we do not fully sample the radial behavior of ct* . We 
have measured a^, within apertures ranging from 2" to 10", and 
the values only vary by < 10% within this range. To reiterate, 
we are not reporting a^, within r^ at this time because we do not 
wish to fold the (large) uncertainties in the structural measure- 
ments into our dispersion measurements. In any case, the radial 
variations in a^, are insignificant on the scales we can currently 
measure. On the other hand, as NGC 2273 demonstrates, two- 
dimensional a^, profiles with higher spatial resolution may well 
reveal a decline in a^, within some galaxy centers. 

Of course, there is the much more serious concern that these 
galaxie s do contain small classical bulges. The recent pa- 
per by iNowak et al.l (12010 ) presents direct evidence that iso- 
lating the classical bulge component in composite systems can 
bring later-type galaxies in line with the local elliptical galaxy 
Mbh - cr* relation. We hope to investigate this concern directly 
in future work. 



7. THE RELATION BETWEEN BH MASS AND BULGE MASS 

Using the subsample of galaxies with SDSS data, we are able 
to investigate the correlation between Mbh and Mbuige- We first 
show the MBH-ibuige relation in Figure 11 a, along with the 
best-fit relation from .Giiltekin et aL (2009). The V-band lu- 
minosities are calculated from the r-band magnitudes and g-r 
colors, using a grid of average galaxy colors derived from the 
Coleman et al. ( 1980) templates. We caution that the sample is 
still small. Furthermore, we are only spatially resolving the in- 
ner ^ 500 pc of these galaxies, while the nuclear spiral in NGC 
2273 is < 300 pc. At the same time, dust extinction is clearly 
significant in all of these systems, and while GALFIT now in- 
cludes dust lanes, the uncertainties are substantial (e.g., B/T 
can vary by nearly a factor of two depending on our treatment). 

With these caveats in mind, we find that the location and scat- 
ter of the megamaser galaxy bulges analyzed so far is consistent 
with that of the early-type inactive sample. We see a slight ten- 
dency for the megamaser galaxies to be overluminous at a fixed 
BH mass. Given that these are spiral galaxies, with bulge col- 
ors that are significantly bluer than typical elliptical galaxies, 
the offset may be attributed to stellar population differences. 
Thus, in Figure 11/? we also show the Mbh -Mbuige relation, 
where the bulge T^ is in ferred from the bu lge g-r color (Table 
2) using the relation of lBelletai] (|2003|) . We find that these 
maser galaxies are consistent with the fiducial Mbh -Mbuige re- 
lation. On the other hand, we find no evide nce that Mbh corre- 
lates with total gala xy luminosity or mass jjahnke et alj|2009t 
lBennertetal.]|2010l §8.2). 
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Fig. 9. — The relation between BH mass and bulge velocity dispersion, including not only the maser sample {open circles) but also the sample from lGiiltekin et al] 
1200911 . Each panel shows fits to a different subsample (see also §6.1). In each panel, literature maser galaxies are highlighted with open grey stars. IC 2560 is 
indicated with a cross and the BH mass eiTor bar is heuristic only. For reference, the A/bh — c* relation of the elliptical galaxy sample (red dashed line) is shown in 
all panels. (Top left): This panel includes all subsamples and all fits, as follows. From Gultekin et al. we show the elliptical (filled red circles), SO (g reen triangles), 
and sp iral galaxies (blue squares). All of the fits are shown as well, including the best fit to the entire sample (thick black solid), the sample of IGiiltekin et alj 
120091 thin black solid), the elliptical (red long-dashed) and non-elliptical (green dot-dashed) subsamples, and to the early (E/SO; magenta long dot-dashed) and 
late-type (blue short dashed) sub-samples. For the purpose of clarity we have omitted upper limits from the figures although they are included in the fits (see §6). 
(Top right): Here we plot the "nonelliptical" sample, consisting of SO (green triangles) and later-type spiral (blue squares) galaxies. The original Giiltekin et al. 
fit to the nonelliptical galaxies (thin green line) and our fit including the maser galaxies (green dot-dashed line) are shown. (Bottom Left): The "late-type" spiral 
galaxies alone (excluding the SO galaxies; blue squares). Best fits to the sample excluding (thin blue solid) and including (blue short-dashed) the maser galaxies are 
shown. (Bottom Right): The "pseudobulge" subsample of the nonelliptical sample is highlighted here with purple asterisks. Note that both SO and later-type spii'al 
galaxies are included, as well as the majority of the maser galaxies. Fits excluding (thin purple solid) and including (purple short-long dashed) the maser galaxies 
are shown. Furthermore, we have performed a second fit (Pseudol) excluding NGC 3245 and NGC 4342, since they have conflicting morphological designations in 
the Hterature. The two galaxies are highlighted with large black triangles, and the alternate fit is shown (dotted purple line). 
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Table 4. Fits 



Sample 


N„ 


Nu 


Q 


13 


e 


<c.o 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 




G09 


49 


18 


8.13±0.07 


4.31±0.38 


0.43±0.05 


^^ 


G09+M 


58 


18 


8.09±0.07 


4.41±0.37 


0.44±0.05 


'>" 2.25 

G 


Elliptical 


25 


2 


8.23±0.08 


4.05±0.40 


0.30±0.05 


NE+M 


32 


16 


7.96±0.13 


3.93 


0.48±0.06 


NE+M 


32 


16 


7.95±0.13 


3.91±0.72 


0.48±0.06 


Early 


38 


6 


8.22±0.07 


3.93±0.37 


0.34±0.05 


b 


Late+M 


19 


12 


7.85±0.29 


3.93 


0.47±0.08 


Late+M 


19 


12 


7.79±0.22 


3.54±1.19 


0.47±0.07 


2 


Pseudo+M 


16 


2 


7.86±0.09 


3.93 


0.29±0.06 




Pseudo+M 


16 


2 


7.88±0.14 


4. 11 ±0.77 


0.29±0.06 




Pseudo2+M 


14 


2 


7.79±0.09 


3.93 


0.26±0.05 




Pseudo2+M 


14 


2 


7.61±0.14 


2.77±0.77 


0.23±0.05 





1.75 



Note. — Col. (1): Fitted Subset. G09=Gultekin et al. 2009; NE= 
Nonelliptical galaxies; Pseudo=Pseudobulges; Pseudo2=Pseudobulge 
sample from G09 excluding NGC 3245 and NGC 4342; M=Megamaser 
galaxies. Col. (2): Number of detections included in the fit. Col. (3): 
Number of upper limits included in the fit. Col. (4): Best-fit intercept. 
Col. (5): Best-fit slope. In cases with no eiTor bars the slope was held 
fixed. Col. (6): Best-fit intrinsic scatter. 

If the Mbh - (y* relation does not hold for these objects, 
while MBH~A/bui g e does , the galaxies ought to be outliers in the 
iFaber & JacksonI (Il976h relation, which we briefly investigate 



7.1. The Faber-Jackson relation 

In principle, examining the location of the maser galaxy 
bulges in the Faber-Jackson relation is an additional tool for 
evaluating bulge morphology. As seen in Figure 12, these five 
galaxie s are cons istent with the Faber-Jackson relation seen in 
the Giiltekin et al. (2009) sample. 

We have argued that a large fraction of the galaxies in our 
sample contain pseudobulges (with the probable exceptions 
of NGC 1194 and NGC 6323). The literatu re on flie Faber- 
Jackson relation in pseudobulges is limited. Korm endyl d 19931) 
find that while the Faber-Jackson relation of elliptical galaxies 
defines an upper envelope for the pseudobulges, they tend to 
scatter toward lo w values of a^ at a fixed B-band luminosity 
(see also Whitmo re et al.|[l979t .Kormendv & Kennicutt 2004). 
The general finding t hat pseudobulg es are more diffuse at a 
given luminosit y (e.g., Caro lldl 19991) i s in th e same sense. On 
the other hand, iGadotti & Kauffmannl (l2009t) rep ort the oppo- 
site, na mely large cr* at a fixed stellar mass. Finallv'Gan da et al] 
(|2009|) report no difference in the //-band Faber-Jackson rela- 
tion of late-type spirals, but they are comparing only to early- 
type spirals, not elliptical galaxies. We await a larger maser 
sample to draw definitive conclusions. 



8. RAMIFICATIONS 

We have demonstrated that there is a population of BHs that 
scatter below the Mbh - cr* relation defined for elliptical galax- 
ies. In this section we explore some of the ramifications of a 
mass-dependent Mbh - cr* relation both for our understanding 
of BH demographics and for our ability to infer BH masses 
from (T* in low-mass galaxies. 
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Fig. 12. — The relation between V— band magnitude and velocity disper- 
sion of the five megamaser galaxies with SDSS imaging {large open circles). 
On their own, the galaxies do not define any notable relation between velocity 
dispersion and magnitude, but they are consistent with the Faber-Jackson rela- 
tion seen in the elliptical and SO galaxies in the lGiiltekin et alj I2009h sample 
(small grey points). 

8.1. The Local Black Hole Mass Function 

It is common practice in the literature to take the observed 
distribution of bulge velocity dispersio ns or magnitudes and 
infer the local BH mass function (e.g.. lYu & Tremaing 120021; 
iMarconi et al.l l2004). If either the scatter or the shape of these 
relations change at low mass, then the shape and amplitude of 
the inferred BH mass function for masses < 10^ M© will also 
change. Increased scatter will tend to broaden the overall dis- 
tribution, thus increasing the relative density of objects at high 
and low mass. The putative offset and change of slope, in con- 
trast, decrease the overall mass density for a fixed distribution 
in cr». At present it is unclear whether the mass functions in- 
ferred from bulge luminosity will change. For this reason, and 
because cr» distributions are not well-measured in the regime of 
interest, we do not present a new BH mass function, but merely 
note that current estimates of the space densities of BHs with 
10'' <Mbh/M0< 10^ are probably overestimates. This is rele- 
vant, for instance, to projected tidal disrup tion rates in present 
and upcoming time-domain surveys (e.g., IStrubbe & Ouataera 
l2002). 



8.2. Normalization of the AGN mass scale 

Currently, the most direct means to estimate BH masses 
in active systems comes from reverberation mapping, in 
which a size scale for the broad-line region is derived from 
the lag between continuum variability and the correspond- 
ing variability from the photoionized broad-line region gas. 
Combining the broad-line region size with the width of an 
emission line yields a "virial" BH mass to within a scal- 
ing factor / that depends on the kinematics and structure 
of the broad-line region. Eventually we hope to measure 
/ in individ ual AGNs and w e are making good progress 
(see, e.g., iBentz et all l2008l 12009c; .Dennev et al., ,20d9l) . 
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Fig. 10. — (a): Large-.scale disk axial ratio versus deviations from Mbh"""* relation; HMs]\=\og{Msi{) — \og[M(at)}. There is no clear correlation between the 
two, suggesting that our cr, measurements are not strongly biased by rotation in the central region, (b): We investigate a possible correlation between the Eddington 
ratio Lbol/^Edd and A Mbh- Again, we see no evidence for a significant correlation between the two, suggesting that deviations from the Mgn — ""» relation are not 
driven by ongoing BH growth. 
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Fig. 1 1 . — (a): The relation between BH mass and bulge luminosity (large open circles) and between BH mass and total galaxy luminosity {small open circles) 
for the megamaser galaxies with SDSS observations. Bulge luminosities are derived by detailed image decomposition of SDSS r— band images using GALFIT (see 
§5, Fig. 7), while th e conversion to V-b and luminosity is calculated as a function of g-r color using the Coleman et al. ( 1980) galaxy templates. The solid line 
shows the relation ofl Giiltekin et al. 1 (2009), while the small grey points are the objects from their sample in order to demonstrate the scatter typical in the inactive 
local samples. The low-mass galaxy anchoring the fit to the galaxies from GUltekin et al. is M32. Note that Gtiltekin et al. fit only E/SO galaxies. Filled black 
circles are AGNs with reverberation-mapping as presented in Bennert et al. (2010) here plotted against total galaxy luminosity, while the squares are local AGNs 
with reverberation mapping that are not included in the Bennert et al. (2010) compilation (§8.2). Finally, the dashed line shows the best-fit relation between Mbh 
and total galaxy luminosity derived by Benne rt et al.. (2010) for AGNs with reverberation-mapping. In these local active samples, including the maser galaxies, total 
galaxy luminosity does not correlate well with BH mass, (b): The relation between BH mass and bulge mass for the same galaxies as in (a). In this case stellar 
masses ar e derived from the r— band magnitude and an estimate of mass-to-Hght ratio in the r— band (T,-) derived from the g-r color using the fitting functions of 
IBell et alj (2003). IC 2560 is indicated with a cross and heuristic BH mass error bars. The upper Hmits (small open circles) are total galaxy masses as derived from 
the r— band magnitude in combination with the T,- for the entire galaxy. Here we show for comparison the fit of Haring & Rix 1 20 04) (solid) and their sample (small 
grey points). Note that while there is substantial overlap between the two inactive samples in (a) and (b), the BH mass measurements are not identical in all cases. 
Again we note that BH mass does not coiTelate with total galaxy mass in the maser galaxies. 
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Fig. 13. — The relation between BH mass and bulge velocity dispersion, 
now including reverberation-mapped objects (grey filled circles). We show 
their virial products (i.e., without any assumption about the broad-line region 
geometry), but we show the true BH masses for the maser galaxies {large open 
circles). IC 2560 is indicated with a cross and the BH mass error is heuristic. 
For reference, the best-fit Mbh~o"* relations using all dynamical masses {thick 
black solid line), just nonelliptical galaxies {thick green dot-dashed line), and 
just late-type galaxies {thick blue dashed line) are shown. We then determine 
the best zero-point offset to match the reverberation-mapped objects with the 
^BH ~ o"* fits to either the nonelliptical galaxies {thin green dot-dashed line) 
or the late-type galaxies {thin blue dashed line). In order to force the active 
galaxies to obey the local Mbh ~ "■* relation for late-type galaxies, one needs 
to multiply the virial pro ducts by a factor / R ^ 3-4, differing from the canoni- 
cal value of / = 5 .5 from lOnkeneTd] f2004l) . 

For the moment, an average / is derived as the BH mass zero- 
point that brings ensembles of low-redshi ft, low-mass AGNs 
into agreeni e nt wi t h the Mbh-ct j, (e..g. iNelson et alj 



2004; 



2008; 



Onken et al.' '2004t iGreene & Hoi t2006b"; 'Shen et al.l 
Woo et al. 2010) or the MsH-ibuige (Ben tz et al. 2009b l) rela 
tions. If the Mbh - cr* relation is intrinsically different for low- 
mass or late-type systems, then / must change as well, since 
the majority of the systems used to derive / have low masses 
and late types dBentz et aljl 2009b). 

In order to illustrate the sense of the shift, we derive / here 
based on our new Mbh - cr* relation for either the nonelliptical 
(including SO galaxies) or the late-type (excluding SO) galax- 
ies. Note that we are ignoring the large uncertainties in the 
measured slope and zeropoint, simply as illustration. Our pro- 
cedure is very simple. We fix the slope of the relation as derived 
from one of the fits above, and then perform a least-squares fit 
for the zeropoint of the active galaxies using the unsealed virial 
products (e.g., v^rJG) as tabulated by Bentz et al. (2009a). 

Starting with the Mbh - cr* relations (Table 4) we derive 
/ = 4.4, 3.4 using the nonelliptical and late-type fit respectively. 
The most recent calibration of /, from Woo et al., is / = 5.2. 
Thus, we infer that / may be a factor of '-^ 1.5 lower than the 
nominal value used in the literature (Fig. 13). In certain con- 
texts a difference of this size matters. Most notably, studies 
looking at potential evolution in the Mbh - cr* or Mbh - Mbuige 
relations wit h cosmic time uniformly adopt th e / value from 
Onken et al. JTreu et al.ll2007t IWoo et al.ll2008l) . However, the 

''in order to calculate a galaxy luminosity w ithout AGN con tamination, we 
V-band magnitude using B-V = 0.8 it 0.2 mag jFukugita eFalliggSi) 



iTreu et al] ( l2007h sample is clearly composed predominantly 
of late- type galaxies. If they were to adopt a lower value of 
/, then the evidence for redshift evolution, which they report 
to be a factor of 2-5 at z = 0.4, could disappear dWoo et al.l 
|2008|) . We note also that the mild offset that is app arent in the 
Mbh ~ o"j, relation of the lo west-mass AGNs (e.g.. lBarth et al.l 
20051 iGreene & Holl2006bl) would also disappear with this re- 
vised scaling. On the other hand, scaling relations at these 
low masses are particularly c omplicated, given t hat we have 
no direc t Mbh measurements dGreene et al.ll2008l) . iBentz et al] 
( 12009 b') find a very well-defined relation between BH mass and 
bulge luminosity for active galaxies with reverberation map- 
ping. With the proper imaging in hand, we will attempt the 
same exercise using the Mbh - Mbuige relation of the maser 
galaxies. 

We wish to make one more point related to active galaxies 
in particular Two recent papers have found that total, rather 
than bulge, luminosity correlates more tightly with Mbh in ac- 
tive systems and that the relation does not evolve with c osmic 
time out to z w 1 dJahnke et al.ll2009l; iBennert et alj|2010h . The 
favored explanation is that secular processes are growing the 
bulge components from the galaxy disks such that the correla- 
tion with total galaxy luminosity is preserved. In Figure lla 
we show the Bennert et al. relation between total galaxy lumi- 
nosity and BH mass measured for local reverberation-mapped 
sources {dashed line). In addition, we include a number of 
nearby reve rberation-mapped galaxies that both Bennert and 
iBentz et al.1 (|2009b) exclude from their samples (NGC 3227, 
NGC 3516, NGC 4051, NGC 4151, NGC 5548, NGC 7469). 

Unlike the Bennert objects, neither the NGC galaxies nor 
the maser galaxies obey a correlation between total galaxy lu- 
minosity and Mbh- Indeed, the lack of correlation with total 
galaxy mass in inactive sampl es was an original motivation t o 
focus on bulge properties (e.g.. lKormendv & Richstorig|1995L , 
Kormendy, J. et al. in preparation). What we do not understand 
at the moment is why, for certain subsets of the reverberation- 
mapped and higher-redshift AGN samples, total galaxy lumi- 
nosity appears to provide a tighter correlation with Mbh- We 
are tempted to blame uncertainties in the zeropoint of the AGN 
masses (see above). However, we do not yet have a full expla- 
nation for this apparent discrepancy. 

8.3. Bias in the Active Sample 

Ideally, we would like to use our results to draw more general 
conclusions about the population of BHs with Mbh~ 10^ Mq 
that live predominantly in spiral galaxies. However, we must 
also consider the possibility that the maser sources represent a 
biased sample. First of all, the selection process may be bi- 
ased by the presence of nuclear activity. Were we dealing with 
distant quasars, we might be selecting on BH mass, while in lo- 
cal inactive systems we select on galaxy properties. However, 
these AGNs are heavily obscured, and at nearly every wave- 
length stellar light dominates the galaxy luminosity on large 
spatial scales. Thus, the masers were selected first by virtue of 
their galaxy luminosity and then based on emission-line ratios 
and we do not expect a large selection bias based on nuclear 
activity. 

A more serious concern is that we do not see a representative 
snapshot of the galaxy population when focusing on accret- 
ing BHs. Active BHs have not necessarily reached their final 

remove the nuclear luminosity as measured by IHo & Peng 120011) and convert to a 



Megamaser Hosts 



17 



masses, and t hus may be gr owing tow ards the i nactive Mbh-ct* 
relation (e.g. jGrupe & Ma thur 2004- iHo et 11 2008; Ki m et all 
1200 8). The BHs presented here are radiating at < 10% of 
their Eddington limits, and thus would need > 1 Gyr to grow 
onto the inactive relation. This is somewhat long compared to 
the ^10^ year lifetimes t ypically estimated for quasars (e.g., 
iMartini & Weinberg|l200ll) . Typical growth times for pseudob- 
ulges , in contrast, are probably many Gyr (e.g., Fish er et all 
|2009|) . Therefore, the notion that these particular BHs knew to 
accrete material at this moment in order to move towards the in- 
active relation seems somewhat contrived. In addition, we plot 
Z\Mbh against Eddington ratio in Figure \Qb. There is no corre- 
lation between these two quantities (Kendall's r = -0.2, prob- 
abihty of no correlation is P = 0.7). We see no evidence that 
the maser galaxies provide a biased tracer of BH-bulge scaling 
relations. 

9. SUMMARY 

We consider the host galaxy properties of a new sample of ac- 
tive galaxies with megamaser disks. Fitting of the maser spots 
yields BH masses with < 15% precision (C. Y. Kuo et al. in 
preparation). BH mass measurements in the maser galaxies 
more than double the sample of galaxies with dynamical BH 
masses of Mbh= 8x1 0^ -7x1 0^ Mq, as well as increasing 
significantly the number of spiral galaxies with dynamical BH 
mass measurements. Our primary results are as follows: 

1 . There is a population of spiral galaxies, probed by the 
maser measurements, with BHs that systematically lie 
below the Mbh - cr* relation of massive elliptical galax- 
ies. As a result, the Mbh-ct, relation fits to later-type 
and lower-mass galaxies display both a larger scatter 
and lower zeropoint than fits to elliptical galaxies alone 
(§6). There is no universal power-law Mbh -cr, relation. 

2. From our limited data set, the galaxies appear to obey 
scalings between Mbh and bulge luminosity or mass. 
On the other hand Mbh does not correlate strongly with 
total galaxy luminosity for the maser galaxies (§7). 

3. To rule out definitively the presence of small classical 
bulges buried in the bars, dust, and young stellar popu- 
lations found at the centers of these galaxies, we require 
spatially resolved kinematics and high-resolution imag- 
ing, which will fully characterize the complicated bulge 
regions of these galaxies. 

4. Changes in the Mbh - cr* relation at low mass translate 
directly into uncertainties in the local BH mass function 
forMBH^ 10^ Mq. Furthermore, they imply ambiguity 
in the overall scaling of BH masses in active galaxies 
based on reverberation mapping. 

We now indulge in speculation about the physical rami- 
fications of our results. Generic arguments about BH self- 
regulation, taken at face value, seem to predi ct a simple 
power-law relationship between Mb h and g^, (e. g., |Silk & ReesI 
[1998; Murray etal. 2005; Hopkins etaD |2006!), although the 
shape and scatt er of the relation may depen d on Mrh (e.g ., 
[Robertson etani2006l; ICroton et al.ir2006: .DeBuhretal.ll20ia) . 
This relation is estabUshed when the BHs grow large enough 
that their accretion power can overcome the galaxy potential. 
On the other hand, the maser galaxies have systematically lower 



BH masses at a fixed cr*. Perhaps because of the quiescent his- 
tory of these spiral galaxies, there has never been a major feed- 
ing episode in the past history of their BHs. If low-mass galax- 
ies do not feed their BHs very effectively, then the BH never 
approaches the "limiting" value set by the ga laxy potential (jHul 
l2008HGreene et al.l2008l; I Younger et al.l2008 ). Presumably the 
same processes that make a galaxy into an elliptical or spiral 
(merging history) also determine the BH feeding efficiency. In 
this scenario, we expect a measurable difference in scatter be- 
tween elliptical and spiral galaxies with Mbh~ 10^ Mq. We 
note that ongoing adaptive-optics-assisted integral-field unit 
observations are exploring the lo w-mass regime for e lliptical 
galaxies with improved statistics (iKrainovic et al.l2009l) . which 
will allow us to perfom this t est. 

Alternatively. iPena (l2007h proposes that BH-bulge scaUng 
relations are tightened at high mass due to multiple minor merg- 
ers, which, by the central limit theore m, yield a tight Mbh ~ cr* 
relation at the present day (see also iHirschmann et al.l 120101; 
iJahnke & Maccioll2010l) . In this picture the only important pa- 
rameter is number of mergers, rather than detailed galaxy mor- 
phology. Here we might expect to find increased scatter in both 
the elliptical and spiral galaxies at Mbh;$ 10^ Mq. It is inter- 
esting to note, in passing, that one might also expect globular 
cluster frequen cy to correlate more tigh tly with Mbh than bulge 
luminosity (Bur kert & Tremainel2010h . Furthermore, we might 
also expect this process to lower MBn/Mbuige over cosmic time 
(e.g.. Greene et al. 2010) . 

This paper represents an initial exploration of scaling rela- 
tions for the maser galaxies, but there are many experiments we 
would like to perform. As noted above, further observations 
are required to determine whether or not there are small classi- 
cal bulges embedded in all of these galaxies, wh ose cr* and M* 
correlate tightly with Mbh (iNowak et alj|2010l) . Integral-field 
spectroscopy would reveal the cr, profile and allow us to deter- 
mine whether there is a transition scale at which we recover a 
tight correlation between Mbh and cr, . Multi-color photometry, 
ideally with HST imaging of the nuclei, will allow us to decom- 
pose the light profile into its constituent parts, measure stellar 
masses for each, and then determine whether the BH mass cor- 
relates best with any particular combination of components. 

There are other potentially interesting BH-bulge scaling re- 
lation projections to investigate as well. For instance, it may 
be that Mbh correlates more tightl y with the dynamical mass 
than cr» or stellar mass (e.g., Haring & Rix 2004; Hop kins et al] 
2007) . There is also a suggestion in the literature tha t Mbh ac- 
tually correlates with the dark matter halo mass (e.g.. lFerraresel 



tually ( 
120021) . 



al though the re are clearly complications in spiral galax- 
ies (e.g.. lHoll2007l) . At present only a couple of the galaxies 
have resolved rotation curves, but we will obtain H I observa- 
tions of the entire sample with the EVLA in order to investigate 
both the gas content and the Mbh-Vc correlation in these galax- 
ies. 

Beyond the present work there are now many suggestions 
that the Mbh - cr* relation may not be the unbroken, low- 
scatte r power-law originally pr oposed (e.g., iGebhardt et al.l 
120001: iFerrarese & Merrittl l2000l) . With the larger samples 
available to date, the scatte r has increa sed by '^ 0.1 dex 
( JTremaine et al. 2002; Gultekin"erai]|2009l) . At the same time, 
new analysis incorporating the im pact of dark matter at large ra- 
dius toebhardt & Thomasl2009 ') has shown that the BH masses 
in the most massive ellipticals may be systematically underes- 
timated by a factor of ^ 2. Finally, new dynamical models in- 
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corporating txiaxiality may well inc rease the scatter in BH mass 
at a fi xed bulge mass as well (e.g.. Ivan den Bosch & de Zeeuwl 
I2OIOI). Given that there is no universal Mbh-ct* relation, it may 
be time to revist the role the BHs play in galaxy evolution. 



Many people have contributed substantively to our thinking 
in the writing of this paper. We first thank the referee for a very 
prompt and thorough report. We gratefully acknowledge use- 
ful conversations with A. Barth, L. Ho, J. Kormendy, G. van de 
Venn, and N. Drory. P. Erwin and J. Mulchaey both provided 
optical images of various galaxies presented here and K. Giil- 
tekin assisted in our implementation of his fitting methodology. 

This research has made use of the NASA/IPAC Extragalac- 
tic Database (NED) which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under contract 
with the National Aeronautics and Space Administration. We 
further acknowledge the usage of the HyperLeda database 
(|http://leda.univ-lyonl .fr). 

Funding for the SDSS and SDSS-Il has been provided by 
the Alfred P. Sloan Foundation, the Participating Institutions, 



the National Science Foundation, the U.S. Department of En- 
ergy, the National Aeronautics and Space Administration, the 
Japanese Monbukagakusho, the Max Planck Society, and the 
Higher Education Funding Council for England. The SDSS 
Web Site is http://www.sdss.org/ 

The SDSS is managed by the Astrophysical Research Con- 
sortium for the Participating Institutions. The Participating In- 
stitutions are the American Museum of Natural History, As- 
trophysical Institute Potsdam, University of Basel, Univer- 
sity of Cambridge, Case Western Reserve University, Uni- 
versity of Chicago, Drexel University, Fermilab, the Insti- 
tute for Advanced Study, the Japan Participation Group, Johns 
Hopkins University, the Joint Institute for Nuclear Astro- 
physics, the Kavli Institute for Particle Astrophysics and Cos- 
mology, the Korean Scientist Group, the Chinese Academy of 
Sciences (LAMOST), Los Alamos National Laboratory, the 
Max-Planck-Institute for Astronomy (MPIA), the Max-Planck- 
Institute for Astrophysics (MPA), New Mexico State Univer- 
sity, Ohio State University, University of Pittsburgh, University 
of Portsmouth, Princeton University, the United States Naval 
Observatory, and the University of Washington. 



200 



150 - 



100 




J I I L_ 



100 



150 

^our (km s-i) 



200 



Fig. 14. — Comparison of stellar velocity dispersion measurements using our internal velocity template stars (cTour) and the Elodie template stars (cEiodie) for the 
IMACS (circles) data. The template stars, while not identical, are matched approximately in spectral type to remove ambiguity caused by template mismatch. We 
measure the differential resolution by fitting the internal template stars with the Elodie templates iMoultaka et al. 2004) and use these values to coiTect the Elodie 
fits for differing instrumental resolution. For demonstration, two additional galaxies are plotted here that are not in the current sample. 



In addition to the B&C data discussed in the primary text, we obtained auxiliary spectra for some of the S outhern targets in one 
Magellan/Baade night. We used the Inamori-Magellan Areal Camera & Spectrograph (IMACS : lDressler et al...2006.) with the 1200 
1 mm"' grating. The primary setting had a tilt of 17.15 with a resulting central wavelength of 4540 A, a spectral range of 3800-5300 
A, and a dispersion of crinsti- ~ 15 km s"'. We also obtained CaT observations for two targets, with a tilt of 32.8, a resulting central 
wavelength of 8340 A, and a spectral range of 7500-9000 A. The observing conditions were not optimal; it was only four days from 
full moon, and the typical seeing ranged from 0."7 to 1 ."4 over the course of the night. 

IMACS consists of an array of four by two CCDs with considerable tilt in the spatial direction. Therefore we used the COSMOS 
code developed by G. Oemler, D. Kelson, and G. Walth to reduce long-slit IMACS data. The COSMOS code has a variety of 
benefits over standard long-slit reductions of the IMACS data. For one thing, because the size of the chip gaps is built into the 
code, we are able to solve for a global (rather than a chip-by-chip) wavelength solution in a straightforward manner. Furthermore, 
the code contains an optical model for the instrument and thus does a better job of modeling curvature in the spatial direction. 
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We perform bias-subtxaction, flat-fielding, dispersion-correction, and rectification within COSMOS and then use standard iraf 
routines for tracing an d extracting one-dimensional spectra, including sky subtraction. Finally, we use IDL routines as described by 
iMatheson et alj (|2008) to perform flux calibration and telluric absorption corrections. 

In addition to the IMACS spectra, three tai-gets ( NGC 1194, NGC 438 8, and NGC 6264) have specti-a in the SDSS database. 
These have a spectral resolution of '^ 70 km s~' (e.g.. lGreene & Holl2006ah and completely cover the G-band region. We use Valdes 
template stars for the SDSS spectra as well as the DIS spectra. 

The velocity dispersion measurements proceeded as above, except that our library of template stars taken with IMA CS is limited 
to G and K giants. Thus we use a large library of template stars provided by the Elodie catalog (Moult aka et al.l l2004). These stars 
were observed with R ^ 42000 or ^ 3 km s"' resolution. In all cases we adopt the following stars as our templates: HD025457 
(F5 V), HD148856 (G8 III), HD185351 (KO III), and HD124897 (K1.5 III). 

We fit our internal velocity dispersion templates with the Elodie stars in order to determine the relative broadening between the two 
systems. In the case of IMACS the spectral resolution varies significantly as a function of wavelength, but in the range 4100 < A < 
4500 A, fits to the IMACS stars with the Elodie templates yield a relative dispersion of 12 ± 0.5 km s"' . As an additional consistency 
check, we take our measured dispersions from the internal velocity stars and compare them to those measured with the Elodie stars 
in the identical wavelength range, where the latter are corrected for the resolution difference in quadrature. The result of that test is 
shown in Figure 14. The agreement is completely satisfactory, with a fractional difference of ((Tour- crEiodie)/cour = -0.02 ± 0.03. 
We are reassured both that there is no significant systematic offset between the two sets of measurements and that the scatter is 
considerably smaller than any other source of uncertainty. 

Overall, we find satisfactory agreement between the IMACS and B&C measurements (Table A5). The one exception is the galaxy 
NGC 3393. In this case, our IMACS data were not taken along the major axis, but rather at PA=12°. Since the IMACS measurement 
is nominally consistent with our minor axis measurement, we suspect this PA difference may at least partially explain the discrepancy. 
Furthermore, we are able to test our flux calibration by comparing the fluxes in [O III] between the two sets of observations. Since 
we can match approximately the slit size and position angle, any discrepancies should predominantly reflect differences in flux 
calibration. It is reassuring that our measurements agree within < 30% in all cases but NGC 2960, where the [O III] line falls on a 
noisy part of the IMACS chip. 
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Note. — Col. (1): Galaxy. Col. (2): Stellar veloc- 
ity dispersion (cr*) within 4" (km s~') along the major 
axis using the region surrounding the G-band. Col. (3): 
Same as (2) but along the minor axis if such data ex- 
ist. Col. (4): u^ within 4" (km s~^) as measured from 
either the SDSS (NGC 1194, NGC 4388, NGC 6264) 
or with our IMACS spectra using the region surround- 
ing the G-band. Col. (5): a^ within ~ 4" from the Ca 
Triplet (km s"' ; 8400-8700A). 
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